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The structure of networks describing interactions be¬ 
tween entities gives significant insights about how these 
systems work. Recently, an approach has been proposed 
to transform a graph into a collection of signals, using 
a multidimensional scaling technique on a distance ma¬ 
trix representing relations between vertices of the graph 
as points in a Euclidean space: coordinates are inter¬ 
preted as components, or signals, indexed by the ver¬ 
tices. In this article, we propose several extensions to 
this approach: We first extend the current methodology, 
enabling us to highlight connections between properties 
of the collection of signals and graph structures, such 
as communities, regularity or randomness, as well as 
combinations of those. A robust inverse transformation 
method is next described, taking into account possible 
changes in the signals compared to original ones. This 
technique uses, in addition to the relationships between 
the points in the Euclidean space, the energy of each sig¬ 
nal, coding the different scales of the graph structure. 
These contributions open up new perspectives by en¬ 
abling processing of graphs through the processing of 
the corresponding collection of signals. A technique of 
denoising of a graph by filtering of the corresponding 
signals is then described, suggesting considerable poten¬ 
tial of the approach. 

1. Introduction 

The development of data processing methods is crucial 
with the recent explosion of data made possible by tech¬ 
nological devices. Many systems, previously inacces¬ 
sible by traditional methods of quantitative analysis by 
lack of information, are now described by huge amounts 
of data, so that the limiting factor is now the absence of 
tools to give them sense. Among these systems, many 
can be represented as networks i.e., a set of relation¬ 


ships between entities. Network theory Il3^ has been 
developed to supply toolboxes, such as for detecting 
communities lITSl . in order to understand the underly¬ 
ing properties of these systems, with many successes. 
More recently, connections between signal processing 
and networks theory have tremendously increased. The 
field of signal processing over networks has been ex¬ 
tensively studied in recent years ll43l with the objective 
to transpose concepts developed in classical signal pro¬ 
cessing, such as Eourier transform or wavelets, in the 
graph domain. These works, based on spectral graph 
theory ifTOl . have led to a growing set of significant re¬ 
sults, among them filtering of graph signals (i.e., signals 
defined over a graph) ll4^ l40l . spectral wavelets ll20ll . 
wavelet filterbanks |[32l l43l |33 [34]| vertex-frequency 
analysis 1441 of graphs signals, multiscale community 
mining using graph wavelets Il46l . or sampling for graph 
signals H |8l |30l [T^ |471 added to considerations about 
uncertainty principle on graphs |[T l[T7ll48ll . 

Eor the study of networks themselves, another ap¬ 
proach linked to signal processing has been considered: 
it consists in a duality between graphs and signals, based 
on the transformation of graphs into signals and recipro¬ 
cally, in order to take advantage of both signal process¬ 
ing and graph theory. With the development of the net¬ 
work science ||33l, some works have proposed to map 
time series into graph objects, aiming at using the wide 
range of measures and tools defined on networks to high¬ 
light properties in the time domain; this approach has 
been successfully used in the analysis of nonlinear time 
series EEa, in the characterization of nonlinear dy¬ 
namics Il5^ l42l [T3]l or in the identification of invariant 
structures l[T2ll . with significant results in applications 
such as heart rate [71 or earthquake [2 analyses. 

Conversely, mapping a graph into time series has been 
less intensively studied. Recently, Weng et al. lISTII pro- 
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posed to explore the structure of scale-free networks us¬ 
ing finite-memory random walks, where the values of 
time series at time t are the degree of the vertex visited 
by the walker at step t. The resulting time series, ob¬ 
tained from different real-world networks, exhibit corre¬ 
lations, linked with the scale-free property of these net¬ 
works. With a similar approach, Campanharo et al. f7] 
proposed a random walk based algorithm to map graphs 
into time series by associating a specific value in fhe 
time domain with vertices, with the particularity that the 
graphs are themselves derived from time series. Girault 
et al. IHtII extended this approach in the case where the 
graph is the object of interest, by using semi-supervised 
learning to map vertices to signal magnitudes such that 
the resulting time series are smooth. Using an alterna¬ 
tive approach, Haraguchi et al. Il26ll and later Shimada et 
al. Bdl proposed a deterministic method based on clas¬ 
sical multidimensional scaling (CMDS) to represent the 
vertices of the graph as a set of points in a Euclidean 
space, where the relations described by the edges are 
represented by distances between points. 

An attractive feature of this latter approach, in com¬ 
parison with random walks, is that it is fully determinis¬ 
tic since for a given graph, its representation in the signal 
domain remains the same. A second interesting point is 
that all information is included in the signals: from them, 
it is possible, under some assumptions, to compute ex¬ 
actly the original graph. Our contribution in this article 
elaborates on preliminary contributions along this line in 
1211, and extends the method proposed in BTI on several 
points: in Sectionthe methodology is examined more 
deeply, and application of the method to a wider variety 
of graph structures is studied through illustrations. A ro¬ 
bust inverse transformation is then proposed in Sectionj^ 
to transform back a collection of signals to a graph, with 
a focus on the case where the collection of signals is de¬ 
graded and is not anymore the direct result of the trans¬ 
formation of a graph. Finally, we develop a compre¬ 
hensive framework in Sectionj^to analyze graphs using 
tools from signal processing, namely spectral analysis, 
with an illustration to the denoising of a graph. Note 
that the present work focuses on static graphs; in previ¬ 
ous works, some applications of the present framework 
have been proposed to study and extract the topology of 
dynamic graphs Il2^l22ll^ . 

Notations Throughout the article, the following no¬ 
tations are adopted. Let Q = (V,f) be a simple undi¬ 
rected and unweighted graph, where V is the set of ver¬ 
tices of size n and S the set of edges of size m. We note 


A its adjacency matrix, whose element aij is equal to 1 
if (i, j) G <S, 0 otherwise, and for i,j G V. The terms 
“signals” or “components” are indistinctly used in the 
following, while avoiding the term time series, to pre¬ 
vent the confusion with the case of dynamic graphs. 


2. From graph to signals 

2.1. Transformation using 
muitidimensionai scaiing 

Shimada et al. BTTI proposed a method to transform a 
graph with n vertices into a collection of signals of n 
points indexed by the vertices of the graph by using mul¬ 
tidimensional scaling (MDS) 0. 

MDS is a set of mathematical techniques used to rep¬ 
resent dissimilarities among pairs of objects as distances 
between points in a multidimensional space whose di¬ 
mension is low. Classical MDS (CMDS) is a partic¬ 
ular case of metric MDS where the dissimilarities are 
assumed to be Euclidean distances. The matrix X of 
coordinates in the low-dimensional space can be com¬ 
puted analytically: Starting with a distance matrix A = 
we first compute a double centering of the 
matrix whose terms are squared: B = — ^ J with 
A(^) = a o a and J = In — where In is the 

identity matrix , a n x n matrix of ones, and o 

the Hadamard product. The CMDS solution is given by 
1 

X = Q+Ajj. with A_|_ a diagonal matrix whose terms 
are the strictly positive eigenvalues of the matrix B 
sorted in an increasing order and Q+ is the matrix of the 
corresponding eigenvectors. The resulting signals are 
the components (or columns) of the matrix X. The y-th 
signal is noted X^^\ An alternative approach to find the 
matrix X consists in solving the following optimization 
problem: minxT(X) = YA^j{5ij - dij{X))‘^ where 
dij{X) is the Euclidean distance between the points i 
and j. An algorithm called SMACOF 0 has been de¬ 
veloped to solve such problem. It is worth to note that 
the solution is not unique, as any rotation, reflection or 
translation of points in the Euclidean space will preserve 
the distances. 

In lim . CMDS is used to transform a graph into sig¬ 
nals by projecting vertices of the graph in a Euclidean 
space, such that distances between these points corre¬ 
spond to relations in the graph. A distance matrix be¬ 
tween the vertices of the graph A describes these rela¬ 
tions: Shimada et al. proposed the following definition 
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from the adjacency matrix A for the distance matrix: 

A = A + w(lnl^-In-A) ( 1 ) 

where w is an arbitrary weight strictly greater than 1. 
This definition does not include a notion of proximity 
between vertices, beyond direct linkage: if the two ver¬ 
tices are connected, the distance is equal to 1, otherwise 
it is equal to w. The remoteness between two vertices, 
which can be measured using for instance the length of 
the shortest path between the two vertices, is not taken 
into account in this definition: two pairs of unlinked ver¬ 
tices will have a distance equal to w, whether they are 
close or not in the graph. One of the advantage of this 
distance matrix, unlike a distance matrix based for in¬ 
stance on the length of the shortest path between the 
vertices, is that it is Euclidean (under some assumptions 
on the value of w, as discussed later in the paper), and 
then there exists an exact solution. Furthermore, the in¬ 
formation about the proximity in the graph is no more 
induced from the graph, but is automatically retrieved 
by the algorithm: if the considered Euclidean space is 
low-dimensional, the distances cannot be well-preserved 
and only the distances representing distant vertices will 
be respected. This enables us to look at, component by 
component, how this approximation is performed, i.e., 
which vertices are roughly considered as neighbors. 

For a given distance matrix, it does not exist necessar¬ 
ily a configuration of points X such that distances be¬ 
tween these points are equal to distances defined in fhe 
disfance mafrix. In fhe case where fhe disfance mafrix 
A is defined by Equafion fhe sfrucfure of fhe graph 
as well as fhe paramefer w has an influence on fhe ex¬ 
istence of a solufion X. This influence has been barely 
sfudied in E6\ and iHTTl . In 1^ . fhe aufhors compare for 
differenl values of w fhe qualify of fhe reconsfrucfion on 
several graphs generated using fhe Waffs-Strogalz model 
and fwo real-world nefworks, by fhe comparison of fheir 
edges sefs. They conclude fhaf fhe value w should re¬ 
side befween 1 (excluded) and 1.01. The same kind 
of approach is followed in BTIl . but restricted to Watts- 
Strogatz model where the probability of rewiring varies 
from 0.01 to 1. Their conclusion is that for n = 400, 
w should be comprised between 1 (excluded) and 1.14, 
and that this upper bound depends on the value of n and 
should be as close to unity as possible, without substan¬ 
tial argument for that. We propose in the following an 
upper bound for the value of w to guarantee that the 
matrix A is Euclidean, i.e., there exists a configuration 
X such that the Euclidean distances between points are 
equal to A. The calculation of this upper bound relies 


on the positive-semidefiniteness of the resulting matrix 
B, which is the covariance matrix of X: B = XX^. 
The study of the eigenvalues of B according to w and 
the adjacency matrix A, detailed in|^ gives an upper 
bound for w: w < with n the number of vertices. 

This result agrees with the partial results obtained in Il26]l 
and m : w should be close to 1, all the greater given 
the number of vertices n. In addition, in practice set¬ 
ting w below this upper bound is not required for most 
graphs, although it ensures that the matrix A will stay 
Euclidean, even with a peculiar topology of the graph. 

Alternative distance matrix could be used instead of 
the one proposed in lldTIl . As described previously, a 
natural measure of closeness between two vertices in a 
graph is the length of the shortest path between the ver¬ 
tices Il45l . Distances also based on similarity between 
adjacent vertex sets are worth considering. These two 
alternative dissimilarities have the advantage to be eas¬ 
ily extended to the case where the graph is weighted i.e., 
each edge has a weight giving its importance. There 
also exist alternative methods to CMDS to transform 
graph into points in a Euclidean space. For instance, the 
method of Eaplacian eigenmaps, proposed by Belkin et 
al. 0, have been used in a context of dimensionality re¬ 
duction and data representation, and is based on a diag- 
onalization of the Eaplacian matrix of the graph. Many 
parallels can be done between the two approaches, even 
if they differ on the final representation of vertices in 
the Euclidean space lOSll . Such a comparison is beyond 
the scope of this article and will not be discussed further 
here. 

2.2. Graph models and theoretical results 

This section introduces models that will be used to gen¬ 
erate illustrations of the method throughout the article, 
as well as theoretical results associated to these models. 


k-ring lattices A k-ring lattice is a graph where each 
vertex i is connected to the vertices — — 

1,..., i — 1, i -I- 1,..., i -I- I — 1, i -I- |}, for k G 
{2,4,..., n}. As Shimada et al. discussed in lldTIl . it 
is immediate to find expected eigenvalues and eigenvec¬ 
tors in this case using circulant matrix theory |[T9ll . 

The computation of eigenvalues of B leads to 
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(a) 2-ring lattice with 60 vertices (2RL-60-10) 




(e) Stochastic block model with 2 communities, 60 vertices, = 
0.7 and= 0.1 (SBM60-2) 
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(g) Erdos-Renyi model with 60 vertices andp = 0.4 (ER60-.4) 




(d) Watts-Strogatz model with 60 vertices, fe = 10 and p — 0.1 

(WS60-10-.1) 
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Component 2 


(f) Stochastic block model with 4 communities, 60 vertices, = 0.9 
and Pi, = 0.01 (SBM60-4) 
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(h) Barbell model with two cliques of 20 vertices and a path between 
them of 20 vertices (BAR60) 


Figure 1: Illustrations on several instances of graph models of the transformation of a graph into a collection of 
signals. All graphs have n = 60 vertices. For each sub-figure, the left plot shows a two-dimensional 
representation of the graph, while the right plot displays the first three components, with high energies, 
and the arbitrarily chosen component 42, a low-energy component. 


4 












































































































\ 9 2t7r 

where a = k + (n — 1 — k)w and C, = e ^ is the 
nth root of the unity. The eigenvectors are given by the 
columns of the Fourier matrix. Details of the calculation 
are given in[B| 

When the eigenvalues are ordered by the value of q, 
the eigenvectors are considered in increasing order of 
frequencies. The components are however sorted ac¬ 
cording to the energy of eigenvalues when applying 
CMDS, which correspond to the sorting according to the 
value of q only when k = 2. In this case, the resulting 
signals are then harmonic oscillations whose frequen¬ 
cies increase as lower-energy components are consid¬ 
ered. When k increases, the components are no longer 
sorted by frequencies. 

Watts-Strogatz model The Watts-Strogatz model 
|[50l has been developed to build graphs with small- 
world property, where the average length of the shortest 
paths between vertices is low compared to the number 
of vertices, while the clustering coefficient is high. This 
property has been highlighted in many systems, such 
as for instance social networks. From a /c-ring lattice, 
each edge (u, v) is rewired with probability p to become 
(u, w), where w is uniformly chosen among other ver¬ 
tices. Shimada et al. give a second-order approxima¬ 
tion of the expected eigenvalues and eigenvectors ac¬ 
cording to the probability p, using perturbation theory. 
They showed that the correlation between the approxi¬ 
mation and the actual signals is high when p is low, and 
decreases when p increases, which is consistent with in¬ 
tuition. 

Erdos-Renyi model A random graph is a graph 
whose set V of vertices is fixed and links befween fhese 
verfices are drawn randomly |[T4ll . Random mafrix fhe- 
ory suggesfs fhaf eigenvectors are random vectors in M”, 
nafurally conjecfured to be disfribufed as i.i.d. Gaus¬ 
sian vecfors im. As for eigenvalues, fhe Wigner resulfs 
ll^ show fhaf fhe maximal eigenvalue is comprised be- 
fween np{l — -^) and np{l -|- -^) while fhe remain¬ 
ing eigenvalues follow a semicircle law on fhe inferval 
giving the energies of signals. 

Stochastic block model A simple stochastic block 
model ||29l is used to generate a graph with communi¬ 
ties. Each of the n vertices is assigned to one of the K 
communities, and edges between each pair of vertices is 
randomly and independently drawn according to prob¬ 
abilities depending on the group of vertices: if the two 


vertices belong to the same community, the probability, 
noted pw, is close to 1 while otherwise, the probability 
between vertices of different groups, noted ph, is lower. 
The settings of p^ and pb lead to customized density of 
edges within and between communities. In lOTII . intu¬ 
itions about the shape of the signals are given, in an ap¬ 
plication of segmentation of images. They suggest that 
the eigenvectors of block matrices are piecewise con¬ 
stant with respect to the communities. 

Barbell model The barbell model liJ| is defined as 
fwo cliques, each confaining ric verfices, linked by a 
pafh wifh np vertices. Based on fhe resulfs defined pre¬ 
viously, fhe shape of fhe eigenvectors is expecfed fo be 
close fo fhe ones in fhe case of sfochaslic block model, 
since each clique is communify, as well as fo fhe ones in 
fhe 2-regular laffices, which is encounfered in fhe pafh 
part of fhe graph. 


2.3. Illustrations 


Figure [T] shows illusfrafions of fhe fransformafion of a 
graph info a collection of signals, on several insfances of 
graph models infroduced in fhe previous secfion, all sef 
wifh n = 60 vertices. For each sub-figure, fhe leff plof 
shows a fwo-dimensional represenfafion of fhe sfudied 
graph, while fhe righf plof displays fhe firsf fhree com- 
ponenfs, wifh high energies, and fhe arbifrarily chosen 
componenf 42, a low-energy componenf. 

In all cases, fhe obfained signals are consisfenf wifh 
fhe expected resulfs described in fhe previous secfion. 


For fhe A:-ring laffices (Figure la for fc = 2 and Fig- 
ure[T^for k = 10), fhe signals are harmonic oscillations, 
grouped in pairs of signals wifh fhe same frequency and 
a difference of phase equal fo For fhe high-energy 
componenfs, fhe value of k does nof have any influence. 
However, fhe frequency of signals differs, as expecfed, 
for lower-energy componenfs, as fhe sorfing of compo¬ 
nenfs is guided by fhe eigenvalues. 

When noise is added, here by considering fhe Waffs- 
Sfrogafz model wifh a probabilify of rewiring equal fo 
p = 0.1 (Figure [Tc] for k = 2 and Figure [7d] for k = 10), 
fhe signals are also noisy, as described by Shimada ef 
al. HTTl . According fo fhe value of k, fhe consequences 
of fhe noise on fhe shape of fhe resulfing signals differ: 
for A: = 10, more edges are rewired, buf fhe sfrucfure 
is barely affecfed since each vertex in linked wifh ifs 10 
nearesf neighbors: The ring is fhen preserved, leading 
fo noisy harmonic oscillations. However, when k = 2, 
one single rewiring causes fhe disappearance of fhe ring: 
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the structure of the graph is then different from a regular 
lattice, leading to high changes in signals. 

Figures and [Tf| show two instances of the stochas¬ 
tic block model, the first one with two communities with 
Pw = 0.7 and ph = 0.1 and the second one with four 
communities and pw = 0.9 and ph = 0.01. An interest¬ 
ing observation is that the high-energy components dis¬ 
play the structure of communities, with noisy plateaus 
corresponding to the dense parts of the graph. It is worth 
noting that the number of relevant high-energy compo¬ 
nents is equal to the number of communities in the graph 
minus one, as one component is sufficient to discrimi¬ 
nate two communities. As for the low-energy compo¬ 
nents, they are noisy signals, corresponding to the struc¬ 
ture inside communities. These random structures cor¬ 
respond to a random graph of type Erdds-Renyi, an in¬ 
stance of which is represented in Figure the result¬ 
ing signals do not exhibit any structure, only looking as 
white noise. 

Finally, Figure [Th| displays an instance of the Barbell 
model, which combines features obtained for the pre¬ 
vious instances: plateaus for the first high-energy com¬ 
ponents, as highlighted for the stochastic block model, 
describing the cliques, harmonic oscillation of the third 
component (only the part which corresponds to a path), 
found for the /c-ring lattice, and finally noisy signal for 
low-energy components, as seen for the Erdds-Renyi 
model. 

These illustrations show the connection between 
graph structure and the resulting signals after transfor¬ 
mation, that will be used in Section [^to study the topol¬ 
ogy of the graph using spectral analysis, and to perform 
standard operations, such as filtering, on graphs. Before 
that, a study of the inverse transformation is proposed, 
as it will enable us to represent a collection of signals in 
the graph domain. 

3. Inverse transformation: From 
signals to graph 

3.1. statement of the problem 

Transformation from signals, or time series, to graph has 
been studied in many applications, as in fT\ or 1351, in 
order to use network theory as a tool for the understand¬ 
ing of time series. These methods are nevertheless of 
no use in our case, because the signals are themselves a 
representation of a specific graph, hence it can be seen 
as a restoration problem. Ignoring this would lead to in¬ 
consistent results between the represented and the recon¬ 


structed graphs. Hence, the inverse transformation shall 
take into account this knowledge to preserve the origi¬ 
nal topology of the underlying graph. By construction 
of the collection of signals X, the perfect retrieval of 
the underlying graph is easily reachable, by considering 
the distances between each point: As built using CMDS, 
these distances represent the adjacency matrix. How¬ 
ever, when X is degraded or modified, for instance by 
filtering (see Section|4^, the distances are no longer di¬ 
rectly the ones computed between vertices, even if they 
stay in the neighborhood of these distances. It is never¬ 
theless not possible to directly retrieve the comprehen¬ 
sive set of links of the graph. This case is yet worth con¬ 
sidering if processing of signals is performed in a goal of 
analysis: altering the signals will also alter the distances, 
preventing the direct inverse transformation. 

We propose in the following section a robust inverse 
transformation, based on the thresholding of distances. 
Two enhancements are discussed to improve the separa¬ 
bility of distances. Finally, we discuss two methods of 
thresholding. 

3.2. Robust inverse transformation 

The two contributions we propose in the following aim 
at addressing the drawback raised in the statement of the 
problem: the two distributions of distances representing 
edges and non-edges, can overlap when the collection of 
signals is degraded, preventing an efficient thresholding. 

Fet us consider a collection of signals X with K com¬ 
ponents. The objective of the inverse transformation 
is to obtain an adjacency matrix A, from the distances 
D{X). This adjacency matrix describes the graph Q = 
(V,f). 

Energy-weighted distances An initial observa¬ 
tion to improve the distinction between distances rep¬ 
resenting edges and those representing non-edges is that 
the energies of components are not taken into account 
in the processing of the distance matrix D{X). The 
high-energy components have though a strong influence 
in the description of the global topology of the graph: 
If the distance between two vertices i and j in a high- 
energy component is high, it means that the two ver¬ 
tices are likely to be distant in the graph. Conversely, 
if the distance in a high-energy component is low, then 
the two vertices are likely to be connected in the graph. 
Neglecting the importance of energies of components 
comes to forget the hierarchy of components in their 
description of the global structure of the graph. One 
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Figure 2: Inverse transformation from a collection of signals obtained after transformation of an instance of the 
stochastic block model with four communities (SBM60-4). The figures show the histograms of distances 
between points iri the Euclidean distances: light blue bars represent distances cor¬ 

responding to edges while dark purple bars represent distances corresponding to non-edges. The dashed 
line represents the threshold obtained using Adapted Otsu’s method (see Algorithm |^, while the dot¬ 
ted line represents the threshold obtained using the number of edges of the original graph. The signals 
are considered either non-degraded or degraded, obtained by removing the lowest-energy component. 
The computation of distances is performed without or with considering energies of components (a) Non- 
degraded case and unweighted distances: the distances are either equal to 1 or w. Discrimination between 
edges and non-edges is trivial, (b) Degraded case and unweighted distances: distribution of distances are 
mixed and no simple thresholding is feasible, (c) Degraded case and weighted distances with a = 2: 
the distributions are separated and thresholding is good (d) Degraded-case and weighted distances with 
a = 4: the distributions are better separated. 


way to add this information is to compute a distance 
weighted by the energies of the components: d{X)ij = 

\lYlk=i fife — XjkY with a > 0, where is the en¬ 
ergy of component k, computed as 
normalized such that Ylik=i ^k ~ parameter a 

controls the importance of the weighting: if a is high, 
the high-energy components have a higher importance 
in the computation of distances compared to the low- 
energy components. Figures|^and|^display the distri¬ 
bution of distances representing edges and those repre¬ 
senting non-edges, using respectively a = 2 and a = 4. 
Compared with the distributions in Figure the two 
distributions are better discriminated when a increases, 
enabling a better retrieval of the underlying graph. The 
right choice of a is nevertheless empirical and depends 
on the collection of signals. 

Sequential update of the adjacency matrix 

Weighted computation of distances enables the inverse 
transformation to take into account the energies of com¬ 
ponents, and produces, as shown in Figure a better 
separation of distances. Energies of components could 
be taken into account in another way by considering 
intermediary adjacency matrices A^, where is the 
adjacency matrix obtained by retaining only the first k 
components using the method described in the introduc¬ 


tion. Thus, these different structures, at different scales, 
are successively considered: the first components de¬ 
scribe the structure at large scale, and adding one-by-one 
the components increase weights for the less dominat¬ 
ing edges. Summing over all these intermediary states 
is based on the assumption that if an edge exists in the 
graph, it will be present in many intermediary states. 
The final adjacency mafrix A will fhen have values be- 
fween 0 and K, and has fo be fhresholded fo obfain a 
binary mafrix. The choice of fhe fhreshold is arbifrar- 
ily sef fo This choice has been guided by empirical 
sfudies, by looking af fhe disfribufion of disfances. Infu- 
ifively, if means fhaf a link has fo be presenf in fwo fhirds 
of fhe adjacency mafrix fo be considered as reliable. The 
algorifhm for fhe sequenfial updafe of fhe adjacency ma¬ 
frix is fhen described in Algorifhm 

Figure gives an example of sequenfial consfrucfion 
of fhe adjacency mafrix A on an insfance of fhe Barbell 
model described in Secfion [23| composed of fwo cliques 
of 20 verfices linked by a pafh of 20 verfices (BAR60). 
Figureshows fhe adjacency mafrix of fhe graph. The 
fransformafion is performed on a degraded collecfion of 
signals, obfained by keeping fhe firsl 10 componenfs of 
fhe original collecfion of signals, over 59 componenfs. 


Figure 3c gives fhe adjacency mafrix A affer fhe sum¬ 


ming of adjacency mafrices obfained by refaining suc¬ 
cessively fhe firsl, fhe firsl fwo and fhe firsl Ihree com- 
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(a) Original adjacency matrix (b) Final adjacency matrix with 

it' = 10 


❖ 


0 10 20 30 40 50 

(c) Adjacency matrix after consider-(d) Adjacency matrix before thresh¬ 
ing the first three components: olding at Step 6 of AlgorithmJ^ 

A = Ai + A.2 + As 



Figure 3: Illustration of the sequential update of the adjacency matrix to retrieve the underlying graph from a 
collection of signals. The original graph is an instance of the Barbell model, composed of two cliques 
of 20 vertices linked between them by a path of 20 vertices (BAR60). The collection of signals X is 
degraded by keeping only the first AT = 10 components. The color codes the intensities, from white to 
black. 


















ALGORITHM 1: Sequential update of the adjacency 
matrix 

Ensure: X the collection of signals 
1 : Define A as an empty adjacency matrix 
2 : for k from 1 to AT do 

3: Compute the adjacency matrix Ak from collec¬ 

tion of signals reduced to the first k components 
4: Add Ak to A 

5: end for 

6 : Set values of A lower than ^ to 0 and greater than 
f tol_ 

7: return A 

ponents: A = Ai + A 2 + A 3 . The three communities 
appear visually, as well as the particularity of the mid¬ 
dle one, which is actually a path. The retrieved topology 
is much more accurate when the sum of adjacency ma¬ 
trices goes on: Figure plots the adjacency matrix at 
the end of the loop, before applying the thresholding in 
Step 6 of Algorithmic A = X]fc=i After the thresh¬ 

olding, the final adjacency mafrix A is binary, as plotted 
in Figure]^ The fopology of fhe obfained graph is very 
close fo fhe original one, even if only a small porfion of 
componenfs is refained. 

ALGORITHM 2: Adapfed Ofsu’s mefhod 
Ensure: d fhe disfribufion of disfances 
1 : for Each value of d sef as fhreshold r do 
2 : Compufe fhe sefs ds and dp as fhe values of d 

respecfively lower and greafer fhan fhe fhreshold 
3: Compute wb and wp as fhe proporfion of ele- 

menfs in respectively dp and dp 
4: Compute fhe variances vb and of respectively 

dp and dp 

5: Store for fhe fhreshold r fhe fofal variance wifhin 

classes tr = wpvp + wpvp. 

6 : end for 

7: Selecf r such fhaf tr is minimum. 

8 : return r 


Selection of the threshold The selection of ihe 
fhreshold is a crucial sfep in fhe differenfiation of dis- 
fances represenfing edges from fhose representing non¬ 
edges. A firsf approach, menfioned in iHdll . is fo pre¬ 
serve fhe number of edges in fhe reconsfrucfed graph: 
fhe selecfion of fhe fhreshold is fhen compufed by con¬ 
sidering fhe smallesf disfances as edges unfil fhe num¬ 


ber of edges in fhe original graph is reached. This ap¬ 
proach could be resfricfive if fhis piece of informafion is 
nof available. More specifically, if fhe degraded collec- 
fion of signals has a high impacf on fhe fopology of fhe 
underlying graph, fhen fhe number of edges in fhe re- 
consfrucfed graph could be significanfly differenf from 
fhe number of edges in fhe original graph. We propose 
in fhe following anofher approach fo fhreshold fhe dis¬ 
fances, analogous fo binarizafion in image processing, 
where a gray-scale image is compressed in a black-and- 
whife image: From a range of gray levels, a fhreshold 
is compufed such fhaf fhe defails of fhe picfure are pre¬ 
served af besf when fhe number of levels is reduced to 2 
(black and while). A well-known algorifhm to perform 
such binarizafion has been proposed by Ofsu |[^ fo seg- 
menl gray-scale picfures: Considering Iwo classes, fhe 
algorifhm finds fhe fhreshold among all possible fhresh- 
olds fhaf minimizes fhe variance wifhin classes. In our 
confexl, fhe disfribufion is composed of all disfances be- 
fween pairs of poinfs, which differs from gray levels, 
since fhe number of possible disfances mighf be equal fo 
fhe number of disfances, while in images, fhe number of 
levels is fixed (for inslance 256 levels for an 8 -bil color 
images, whafever fhe size of fhe image). We propose an 
adapfalion of fhe Ofsu’s mefhod in Algorifhm 

The choice of fhe number of edges in fhe consfrucled 
graph is fhen guided by fhe obfained disfances, and is nof 
dependenl of fhe original slrucfure of fhe graph. Figure]^ 
shows for fhe four histograms of disfances fhe obfained 
fhresholding using fhe number of edges (dotted line) and 
fhe adapfed Olsu’s mefhod (dashed line). 

3.3. Performance 

Experimental setup The performance of fhe robusf 
inverse fransformafion using fhe proposed enhancemenfs 
is evaluafed in fhis section. A fhorough evaluafion is al- 
mosf impossible: if fhe collecfion of signals is direcfly 
obfained using a fransformafion from graphs fo signals 
as infroduced in Section fhen fhe inverse fransforma¬ 
fion is immediate and exacf, and if does nof require any 
sophisficafed mefhod. Conversely, if fhe collecfion of 
signals is nof fhe represenfafion of a graph, or if fhe sig¬ 
nals are modified, fhere is no indicafion of fhe correcf 
graph which is described af besf by fhis collecfion. We 
propose nonefheless fo assess fhe mefhod by compari¬ 
son of a graph wifh reconsfrucfed versions affer applying 
perturbations on fhe collection of signals. 

Lef us consider a graph Q = (V,<?), described by 
a collecfion of signals X. The signals are perturbed 
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Quality of reconstruction Quality of reconstruction Quality of reconstruction Quality 



(a) RL60-2 (b) RL60-10 



(c) WS60-2-.1 (d) WS60-2-.1 



(e) SBM60-2 


(f) SBM60-4 



(g) ER60-.4 (h) BAR60 


Figure 4: Results of the inverse transformation of a degraded collection of signals, obtained by retaining only a 
portion of the components. For each number of retained components, a quality of reconstruction is 
measured by comparing the obtained graph with the original one. Each sub-figure shows the results for 
one instance, whose name refers to those given in Figure [T] For each sub-figure, the left plot shows the 
results using thresholding based on the number of edges, while the right plot shows the results using 
the Adapted Otsu’s method. Four configurations of the inverse transformation are studied, none: no 
enhancement is performed; W: weighted computation of distances only with a = 4; S sequential update 
of the adjacency matrix only; WS both enhancements. 
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Quality of reconstruction Quality of reconstruction Quality of reconstruction Quality 



(a) RL60-2 (b) RL60-10 




(c) WS60-2-.1 


(d) WS60-2-.1 



(e) SBM60-2 


(f) SBM60-4 



(g) ER60-.4 (h) BAR60 


Figure 5: Results of the inverse transformation of a degraded collection of signals, obtained by adding Gaussian 
noise with mean 0 and variance ranged from 0% to 50% of the maximal value of the signals. For each 
value of a, a quality of reconstruction is measured by comparing the obtained graph with the original one. 
Each sub-figure shows the results for one instance, whose name refers to those given in Figure [T] For 
each sub-figure, fhe leff plof shows fhe resulfs using fhresholding based on fhe number of edges, while 
fhe righf plof shows fhe resulfs using fhe Adapfed Ofsu’s mefhod. Four configurations of fhe inverse 
fransformafion are sfudied. none: no enhancemenf is performed; W: weighfed compufafion of disfances 
only wifh a = 4; S sequenfial updafe of fhe adjacency mafrix only; WS bofh enhancemenfs. 
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to obtain a degraded collection of signals X, whose 
inverse transformation gives a graph Q = (V,£^). G 
and G are compared by using an index of similarity, 
noted Q{G,G), based on the Jaccard index ll27ll defined 
as Q = , which compares the ratio of common 

edges between G and G over the total number of edges. 
Q{G,G) is comprised between 0, if all the edges are dif¬ 
ferent, and 1, if both graphs have the same edge set. 

Two perturbations of the matrix X are studied. A 
first perturbation consists of removing low-energy com¬ 
ponents of a collection of signals: Starting from the first 
component, the components are successively added one 
by one in decreasing order of energy, until the collection 
is complete. The second perturbation consists in adding 
noise to the components: a Gaussian noise is added with 
mean 0 and variance a, ranged from 0% to 50% of the 
maximal value of the signals. 

The inverse transformation is performed using the dif¬ 
ferent enhancements of the method, in order to assess 
the contributions of each one. Four cases are then dis¬ 
tinguished, according to the presence or absence of the 
weights in the computation of distances (with a = 4), 
and of the sequential update of the adjacency matrix: 
none (no enhancement), W (weighted computation of 
distances only), S (sequential update of the adjacency 
matrix only) and WS (both enhancements). Besides, the 
two methods of thresholding are tested separately in sim¬ 
ilar conditions. The experiments are performed on in¬ 
stances of graph models described in Section]^ in order 
to browse a wide variety of graph structures. 

Results Figures]^ andshow the results of the eval¬ 
uation of the inverse transformation. Each sub-figure 
shows the results for one instance, whose name refers 
to those given in Figure [T] Since the Adapted Otsu’s 
method is not itself an enhancement but an alterna¬ 
tive method to find a threshold of distances, the results 
have been divided for readability in two plots: For each 
sub-figure, the left plot shows the results without using 
thresholding based on the number of edges, while the 
right plot shows the results using the Adapted Otsu’s 
method. 

Let us consider the case where the perturbation of the 
collection of signals X is obtained by retaining only a 
reduced number of components (Figure]^. We first fo¬ 
cus on the left plots of each sub-figure, where the thresh¬ 
old is chosen by using the number of edges. Two con¬ 
figurations emerge from these results: W and WS. The 
gain in quality of reconstruction when using weighted 


distances is obvious, and reveals that this feature is the 
best contribution in the retrieval of the original graph. 
Adding S in combination with W does not affect sig¬ 
nificantly the score obtained with W alone, but with¬ 
out any improvement of the reconstruction, except for 
the instance (BAR60) in Figure In this case, only 


the combination of both enhancements gives good re¬ 
sults. The comparison of the results in left and right 
plots shows that the choice of method of thresholding 
has a slight impact on the results when using the inverse 
transformation in configuration WS: the quality of re¬ 
construction is quite similar in both cases. The Adapted 
Otsu’s method seems however less efficient when the 
number of retained components is low. An immediate 
explanation is that considering the number of edges of 
the original graph in the reconstruction gives undoubt¬ 
edly a higher similarity between the reconstructed and 
the original graph than retaining only distances based on 
a discrimination, as it is performed in the Adapted Otsu’s 
method. Considering nonetheless that the latter method 
gives bad results when the number of retained compo¬ 
nents is low might be a cursory glance: with only few 
retained components, the actual graph described by X 
is likely to be quite far from the one described by X. 
The Adapted Otsu’s method gives then a graph G which 
is not similar to G, but which reflects the actual number 
of vertices. 


Figure]^ shows the results for the second experiment, 
where the degraded collection of signals X is obtained 
by adding Gaussian noise to X. Unexpectedly, the bal¬ 
ance between the importance of W and S is the reverse: 
using enhancements S alone leads to better results than 
using W alone. The combination of both enhancements 
is nevertheless still the best method to reconstruct the 
graph from degraded collection of signals. The compari¬ 
son between the two methods of thresholding is the same 
as in the previous case: Generally, results are a bit better 
when using the number of edges to select the threshold, 
especially when the perturbation is high. 

To conclude this section about the inverse transforma¬ 
tion, two contributions have been done to the method of 
inverse transformation to highly improve the quality of 
restoration from a degraded collection of signals. These 
results also show that not having the knowledge about 
the original graph does not penalize the inverse transfor¬ 
mation, as it is possible to retrieve the correct number of 
edges from the signals. Besides, results give reasons to 
believe that the Adapted Otsu’s method leads to a graph 
which is closer to the actual graph represented by the 
perturbed collection of signals. 
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4. Analysis of signals representing 
graphs 

The last section of this article is dedicated to the anal¬ 
ysis of the collection of signals obtained after transfor¬ 
mation from a graph. The aim of this part is to show 
how spectral analysis can be used to identify specific 
graph structures. As seen in Figure [T] signals present 
specific shapes which can be linked wifh fhe graph sfruc- 
fure. Characterizing fhese shapes using specfral analy¬ 
sis enables us fo associafe frequency pafferns wifh graph 
topology. Before fhaf, we infroduce fhe mefhod of anal¬ 
ysis by addressing an issue relafed fo fhe indexation of 
signals. 


der fo find a collection of independenf pafhs wifh respecf 
fo fhe local sfrucfure of fhe graph, while fhe second step 
defermines fhe besf way fo arrange fhe pafhs such fhaf 
fhe objecfive funcfion of fhe cyclic bandwidfh sum prob¬ 
lem is minimized. Defails of fhe algorifhm and resulfs 
abouf fhe consisfency befween fhe obfained vertex or¬ 
dering and fhe fopology of fhe graph are covered in ESll . 

Figurej^shows fwo examples of Iransformalion of fwo 
parficular nefworks into signals. The resulfing collection 
of signals is indexed by fhe verfices. Three signals are 
displayed on each subplof as follows: Random ordering 
of vertices (Lefl), suifable ordering of verfices (Righf). 
If clearly shows fhaf fhe ordering of fhe vertices has an 
influence on fhe specfral properties of fhe signals. 


4.1. Indexation of signals 

So far, fhe imporfance of indexafion of signals has been 
hidden since only relafions, described by fhe Euclidean 
disfances, have been of inferesf. Hence, fhe order in 
which we considered fhe verfices in fhe inverse frans- 
formafion has no significance. However, fhis order is 
essenfial fo sfudy some specfs of fhe signals, especially 
when using specfral analysis of fhe signals, as fhey are 
indexed by fhis verfex ordering. In Figure [T] fhe verfex 
ordering has been suifably defined fo highlighf specific 
shapes of signals: looking af fhe numbering of vertices in 
fhe graph represenfafion shows fhaf numbers closely fol¬ 
low fhe fopology of fhe graph, in a “nafural” order. Or¬ 
dering randomly fhe verfices does nof change fhe value 
assigned fo each vertex, buf would lead fo abrupf varia¬ 
tions in fhe represenfafion of signals: Specific frequency 
properfies, clearly observable in signals, will no longer 
be visible. Unforlunalely, fhe suifable ordering is usu¬ 
ally nof available, especially when dealing wifh real- 
world graphs. To address fhis issue, we proposed in 
Il25]l fo find a verfex ordering fhaf refiecfs fhe fopology 
of fhe underlying graph, based on fhe following assump¬ 
tion: if fwo verfices are close in fhe graph (by consider¬ 
ing for insfance fhe lengfh of fhe shortesf pafh befween 
fhem), fhey have fo be also close in fhe ordering. The 
mefhod consisfs of fhe sfudy of a relafed labeling prob¬ 
lem, called cyclic bandwidfh sum problem Il28ll . defined 
in fhe more general framework of graph labeling prob¬ 
lems f9]. These problems seek a mapping from vertices 
to infegers, in such a way fhaf an objecfive funcfion, de¬ 
fined as fhe sum of fhe disfances in fhe verfex ordering 
befween all pairs of connected verfices, is minimized. 
The heuristic we proposed in |[25l consisfs of a fwo-sfep 
algorifhm. The firsf step performs local searches in or- 


4.2. Connection between frequency 
patterns of signals and graph 
structures 


Specfral analysis is performed using sfandard signal pro¬ 
cessing mefhods: Eel a collecfion X oi K signals in¬ 
dexed by n vertices, fhe specfrum S gives fhe complex 
Eourier coefficienls whose elemenfs are obfained by ap¬ 
plying fhe Eourier Iransform on each of fhe K compo- 
nenfs of X: s^f = estimated, for positive 

frequencies, on F = § -h 1 bins, F being fhe Eourier 
iransform and k G {l,...,iT}. Erom fhe specfrum S, 
fhe magnifudes M of each frequency / for each com- 
ponenf read as: M{k,f) = |sfc/|. The mafrix M is 
sfudied as a frequency-componenf map, exhibiting paf- 
fems in direcf relafion wifh fhe fopology of fhe underly¬ 
ing graph. The phases of signals are stored in a mafrix 
cf) fo be used in fhe inverse Eourier fransformafion, when 
fhe collecfion of signals has to be refrieved from M. 

Eig. shows fhe frequency pafferns obfained for fhe 
graphs defined in Eigure[T] Each graph highlighfs a spe¬ 
cific frequency paffern, linked wifh ifs topology. Regular 
fe-laffices display single-frequency componenfs, whose 
order depends on fhe value of k: when k is higher fhan 
2 , fhe sorting of componenfs is no longer consisfenf wifh 
fhe increasing order of frequencies, as described in Sec¬ 
tion 2.2 (Eigures |7^ and 7b I. Adding noise to fhe 
graph also affecls fhe pafferns, as can be observed in 
Eigures|7^and|^ Graphs wifh communities are associ- 
afed wifh highly-localized fealures: high-energy compo- 
nenfs, formed by plafeaus corresponding fo fhe commu¬ 
nities, appear wifh high magnifudes for low frequencies 
in fhe firsf componenfs (Eigures]^ and [^. These latter 
patterns are visible in fhe frequency-componenf map of 
fhe Barbell model in Eigurej^ denofing fhe combina- 
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Frequency Frequency 



(a) 2-ring lattice 



Vertices Vertices 


(b) Stochastic block model with 4 communities 


Figure 6: Examples of transformation of two particular networks into signals, indexed by the vertices. The resulting 
collection of signals is indexed by the vertices. Three signals are displayed on each subplot as follows: 
Random ordering of vertices (Left), suitable ordering of vertices (Right). 




(a) RL60-2 (b) RL60-10 (c) WS60-2-.1 (d) WS60-10-.1 



0 10 20 30 40 50 0 10 20 30 40 50 0 10 20 30 40 50 0 10 20 30 40 50 

Components Components Components Components 


(e) SBM60-2 (f) SBM60-4 (g) ER60-.4 (h) BAR60 

Figure 7: Frequency patterns obtained on instances of graphs defined in Section]^ Each pattern can be linked with 
the topology of the underlying graph. The color codes the intensities, from white to black. 
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Figure 8: This diagram describes the process implemented to perform filtering of graphs using duality between 
graph and signals. The boxes describe the objects, while arrows represent the operations. The procedure 
is the following: From a graph Q, a collection of signals X is obtained by performing a good indexation 
of vertices and applying the CMDS method. A spectral analysis of the signals is then performed, and 
a matrix of magnitudes M. is obtained. In the frequency domain, frequencies are truncated, leading 
to a matrix Af, which gives a new collection of signals Y by inverse Fourier transform. This process 
comes to apply a crude low-pass filtering on X. A graph 'K is then computed from X, using the inverse 
transformation defined in Sectionj^ Each object belongs to a specific domain, either graph domain, time 
domain or frequency domain. Operations performs some processing of objects, such as filtering, as well 
as the mapping between the three domains. 


tion of both structures in the graph. Finally, the mag¬ 
nitudes of random graphs, in Figure]^ do not look like 
white noise spectra as the signals are re-indexed, and this 
removes the i.i.d. property and adds higher magnitudes 
on low frequencies for the first components than ex¬ 
pected. Nonetheless, there is no specific pattern, which 
is consistent with the random structure of this graph. 

4.3. Application to filtering of graphs 

Throughout this article, we described tools to transform 
a graph into a collection of signals, analyze these signals 
to describe the graph, and finally transform back these 
signals into a graph, even if the collection of signals is 
not exactly the same. We propose here an application 
of these tools to perform filtering of graphs: in the same 
way as classical signals, graphs are often built from mea¬ 
surements, inducing noise, i.e., presence of undesirable 
edges or absence of existing edges. Figure proposes 
a framework to use duality between graph and signals 
as a way to easily perform filtering on graph by filtering 
signals representing the graph. 

The procedure is the following: From the matrix of 
magnitudes iW obtained from a graph Q, a matrix N 
with the same shape is derived by retaining only the 
lowest frequencies. Using the matrix of phases a 
degraded collection of signals Y is obtained by inverse 


Fourier transform. This process is tantamount to apply¬ 
ing a crude low-pass filter on X. A graph Ti is then 
computed from Y, using the weighted distances with 
a = 2, the sequential update of the adjacency matrix 
and a thresholding using the Adapted Otsu’s method, as 
described in Section |3l 

Figurej^shows two examples of graph denoising, with 
on the left the adjacency matrix of the graph Q, and 
on the right the adjacency matrix of the graph T-L, ob¬ 
tained after applying the process described in Figure 
The first example is an instance of the Watts-Strogatz 
model (Figure |^, corresponding to a noisy fc-ring lat¬ 
tice (WS60-10-.1). The five lowest frequencies have 
been retained, and this leads to strengthening the diago¬ 
nal, characterizing A:-ring lattices. Besides, off-diagonal 
edges have been removed. The second example (Fig¬ 
ure |9^ is an instance of the stochastic block model with 
communities (SBM60-2). This can be viewed in this 
case as two cliques whose random edges inside com¬ 
munities have been removed while random edges be¬ 
tween communities have been added. After denoising, 
where the ten lowest frequencies of the matrix M have 
been retained, the missing edges inside communities 
are retrieved, while edges between them are removed: 
the communities are better-defined, appearing as actual 
cliques. These two examples highlight that denoising 
signals representing the graph is tantamount to denois- 
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Figure 9: Illustrations of denoising of graphs, by applying the process described in Figure]^ The left plot gives the 
adjacency matrix of the graph Q, before filtering, while the right plot gives the adjacency matrix of the 
graph after filtering, (a) The five lowesf frequencies have been refained (b) The fen lowesf frequencies 
have been refained. 


ing fhe graph ifself. 


and only if B is posifive definife: 


5. Conclusion 

In fhis arficle, we proposed a framework fo sfudy graphs 
in fhe classical domain of signals, enabling fhe analysis 
fo lake advanlage of basic signal processing fools. From 
fhe melhod of Iransformafion from graph fo signals pro¬ 
posed by Shimada ef al. iHTI . we discussed exfensions 
fo sfudy fhe oblained signals and indirecfly, fhe graph 
ilself, using a robusl inverse Iransformafion. Denoising 
of graphs has been proposed fo evaluale fhe melhod as a 
firsl pofenfial applicafion. The resulls we oblained sug- 
gesl lhaf fhe framework developed provides a connec¬ 
tion belween signal processing operalions and modifica¬ 
tions of graphs. This paves fhe road fo new approach for 
fhe analysis of real-world nelworks, using more sophis- 
licafed signal processing fools, such as adaptive tillering. 
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(z, Bz) > 0 for all vectors z G M” (3) 

or equivalenfly, if and only if is condilionally neg- 
alive detinile: 

n 

(z, A(2 )z) <0 Vz G M” such lhaf ^ = 0 (4) 

From fhe definition of A in Eq.[T] we have: 

(z,A(2);z) = (2,Az) (5) 

-f Inln^) - {Z, InZ) - {z, Az)) 

= (z, Az) - m^((z, z) -f (z, Az)) 


Two cases can be Ihen distinguished: 

1. If (z, Az) > —(z, z), Ihen 

_ 

(z, Az) -f (z,z) 
which is hold as tt; > 1 . 


( 6 ) 
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6 of fhe region Rhone-Alpes. 


A. Upper bound for the parameter 

w 


2. If (z, Az) < —{z, z), Ihen 


< 


{z,Az) 

{z,Az) + (z,z) 


(V) 


The upper bound depends on fhe adjacency malrix A, 
i.e., on fhe sfruclure of fhe graph. To have an idea of 
a suifable value of w, lei us define A and z such lhaf 
(z, Az) is minimal. A is defined as fhe adjacency ma- 
frix of a graph wilh n vertices, wilh n even, such lhaf 


We propose here a ralionale for fhe choice of w, based 
on fhe sfudy of malrix eigenvalues: As mentioned in fhe 
work of Gower |[T8]| . A is exaclly relrieved from X if 
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Qij = 1 if and only if i and j do not belong in the same 
subset among {1,..., and {| + 1,..., n}. A is then 
a 4-block matrix, with the bottom-left block and the top- 
right block equal to 1. As for z, it is equal to —1 for the 

first half of the vector and 1 for the last half of the vector: 

2 

z = [—1, —1,..., 1,1]. (z, Az) is then equal to 
while (z,z) = n. Hence, we obtain an approximation 


for the upper bound of w\ w < 



B. Theoretical results for the A:-ring 
lattices 


A A:-ring lattice is a graph where each vertex i is con¬ 
nected to the vertices — — — 

1,..., z -I- I — 1, f -I- |}, for fc G {2, 4, .. ., n}. As 
Shimada et al. discussed in ED, it is immediate to 
find expecfed eigenvalues and eigenvecfors in fhis case 
using circulanf mafrix fheory lIT^ . We explicif addi¬ 
tionally here fhe connecfion befween fhe paramefers w 
and k and fhe resulting signals. Any circulanf ma¬ 
frix C has ifs eigenvalues A given \/q G {0,n — 1} 
by Ag = where c is fhe circulanf vec¬ 
tor of C and ^ is fhe nfh roof of fhe unify. 

As for eigenvecfors, fhey are given \/q G [0, n — 1] 
by Vq = ^/n[l, .. ., corresponding to 

the columns of the Fourier matrix noted These 

eigenvalues and eigenvectors appear as complex con¬ 
jugate pairs, namely Ag = Xn-q and Vg = Vn-q for 
g 7 ^ 0. As we consider symmetric matrices, the eigen¬ 
values are real and double (Ag = Xn-q for q > 0 ) and 
the corresponding eigenvectors are the real and imag¬ 
inary parts of Vq, normalized by y/2 to obtain an or¬ 
thonormal matrix: Vq = = \/ 2 cos(^) 

and Vn-q = \/29(i^^'^^) = \/ 2 sin(^), corresponding 
to harmonic oscillations. If n is even, A| is single and 


the corresponding eigenvector is not normalized by 1 / 2 . 

Starting from the circulant vector S with three 
values 0, 1 and w, the circulant vector b is defined 
by bi = —^[Sf — ^], wifh a = k + {n — 1 — k)w‘^. 
The vector b is fhen complefely defined by fhree 
values: ^ when i = 0 , -f(l - ^) when 

z G {1,... U {n - I,..., n - 1} and 

' k 


— >jg _ 

— 2n Z^j=0 


when z G {| -|- 1,..., n — | — 1}. The compulation 
of eigenvalues of B leads fo Ag 

5 (ea. c”+E”:A| c=’ +E”;!;,' c* 

From fhese eigenvalues, fhe obfained signals are similar 
to fhose obfained by slandard melhods of diagonal- 


izalion of malrices, up fo rolalion, reflection and 
franslafion. 
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